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Abstract. The recent detection of TeV gamma-radiation from the direction of 
the Galactic Center within several arc-minutes around Sgr A* is the first model- 
independent evidence of existence of high energy particle accelerator(s) in the central 
10 pc region of our Galaxy. This is an extraordinary site that harbours many remark¬ 
able objects with the compact radio source Sgr A* - a hypothetical super-massive 
black hole (SMBH) - in the dynamical center of the Galaxy. Here we explore the 
possible direct and indirect links of the reported TeV emission to the SMBH. We 
show that at least three y-ray production scenarios that take place close to the 
event horizon of the SMBH can explain the reported TeV fluxes. An alternative (or 
additional) channel of TeV radiation is related to the run-away protons accelerated 
in Sgr A*. Quasi-continuous injection of relativistic protons into the surrounding 
dense gas environment initiates detectable high energy gamma-ray emission. The 
absolute flux and the energy spectrum of this radiation component strongly depend 
on the history of particle injection and the character of diffusion of protons during 
the last 10® years. For a reasonable combination of a few model parameters, one can 
explain the detected gamma-ray flux solely by this diffuse component. 
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1. Introduction 


The recent detection of TeV gamma-ray emission from the direction 
of the Galactic Center by three independent groups, CANGAROO 
(Tsuchiya et ah, 2004), Whipple (Kosack et ah, 2004) and HESS (Aha¬ 
ronian et ah, 2004), is a remarkable result which will have a strong 
impact on our understanding of high energy processes in the central 
region of the Galaxy. The localisation of the TeV signal by HESS 
within a few angular minutes indicates that the gamma-ray source(s) 
is (are) located in the central < 10 pc region. Among the possible sites 
of production of TeV gamma-rays are Sgr A*, the young supernova 
remnant Sgr A East, the Dark Matter Halo, and finally the whole 
diffuse 10 pc region. It is possible that some of these potential gamma- 
ray production sites comparably contribute to the observed TeV flux. 
Moreover, the same source could be responsible for two different compo¬ 
nents of radiation: (i) the direct gamma-ray component produced inside 
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the particle accelerator, and (ii) the indirect component produced by 
runaway protons which are accelerated in the same source, but later 
injected into the surrounding dense gas environment. 

Below we discuss these two components of TeV radiation in the 
context of acceleration of particles in the proximity of the SMBH and 
their radiation both inside and outside of Sgr A*. 


2. Production of high energy gamma-rays in Sgr A* 

The temporal and spectral features of Sgr A* are unusual and, as a 
whole, different from other compact galactic and extragalactic sources 
containing black holes. This concerns, first of all, the extraordinary low 
luminosity of Sgr A*. In addition the other important astrophysical 
implications, the low luminosity of Sgr A* has a dramatic effect on the 
visibility of the source in gamma-rays (Aharonian and Neronov, 2005). 

2.1. Transparency of Sgr A* for high energy gamma-rays 

Because of internal photon-photon pair production, the high energy 
gamma-ray emission of these objects (both of stellar mass and super- 
massive BHs) is generally suppressed, and consequently the unique 
information on possible particle acceleration processes near the event 
horizon of the BH is essentially lost. But this is not the case of the super- 
massive BH located at the dynamical center of our Galaxy (Sgr A*), 
which thanks to its extraordinary low bolometric luminosity (< 10 “®LEdd) 
is transparent for very high energy gamma-rays. It is seen from Fig. 1 
that indeed up to 10 TeV the source is transparent for gamma-rays even 
if one assumes that gamma-rays are produced within 2Rg. Note that the 
decrease of the 77 ^ e'’'e“ cross-section well above the pair production 
threshold makes the source again transparent at E r\j eV 

High energy gamma-rays from compact regions close to the event 
horizon of the SMBH can be produced in various ways due to acceler¬ 
ation of protons and/or electrons and their interactions with ambient 
magnetic and radiation fields, as well as with the thermal plasma. 

2.2. Synghrotron and gurvature radiation of protons 

Synchrotron and curvature radiation of protons are important pro¬ 
cesses in the so-called extreme accelerators (Aharonian et ah, 2002) 
where particles are accelerated at the maximum possible rate, E = eB. 
However, even for an “ideal” combination of parameters allowing the 
most favourable acceleration/cooling regime, the characteristic energy 
of synchrotron radiation is limited by Cmax = (9/4)Q;/^mpC^ ~ 0.3 TeV 
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Figure 1. Attenuation of gamma-rays in Sgr A* due to internal photon-photon pair 
production dominated by interactions of high energy gamma-rays with radiation of 
the compact infrared source. Two solid curves marked “ 7 ” are calculated assuming 
that the infrared emission of Sgr A* is produced within lOiZg and 2Rg around the 
central black hole of mass 3 x 10® M©. The curve marked “n” shows attenuation of 
the neutron flux exp (—d/A), where A « 10(15/10^® eV) kpc is the decay mean free 
path of a neutron of energy E, and d = 8 kpc is the distance to the Galactic Center. 


(Aharonian, 2000). This implies that the proton-synchrotron radiation 
cannot explain the flux observed from the direction of GC up to several 
TeV, unless the radiation takes place in a source moving towards the 
observer with bulk motion Lorentz factor T > 10. Another possibility 
could be if the proton-acceleration and synchrotron gamma-ray pro¬ 
duction regions are separated, e.g. when protons are accelerated in a 
regular B-field while moving along field lines, and later are injected 
into a region of chaotic magnetic field. 

In this scenario we should expect 7 -rays also from the regular B- 
field region - due to the proton curvature radiation. The contribution 
of latter in the high energy radiation of SMBHs could be quite signif¬ 
icant (Levinson, 2000). The curvature radiation of protons can extend 
to Cmax = 3Ep/2m^R ~ 0.2(5/10^ )^/^TeV (hereafter all estimates 
correspond to the black-hole mass 3 x 10® Mq). Formally, this equation 
allows extension of the spectrum to 10 TeV, provided that the magnetic 
field exceeds B ~ 10® G. However, such a strong field would make the 
source opaque for TeV gamma rays (Aharonian and Neronov, 2005). 
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2.2.1. Photo-meson interactions 

The protons accelerated in the region close to the event horizon of 
SMBH to energies E ~ 10^® eV, start to interact with soft photons 
of the compact infrared source located at ~ 10i?g. Despite the low 
luminosity of the source, the density of infrared photons appears suffi¬ 
ciently high for reasonably effective photo-meson interactions. Indeed, 
the mean free path of protons through the photon field is estimated 

~ ~ cm)^ cm. This means that ap¬ 

proximately R/Ap-^ ~ 0.01 fraction of the energy of protons is converted 
into secondary particles (neutrinos, photons, electrons). 

While neutrinos and neutrons, as well as gamma-rays of energies 
below 10^^ eV escape freely the emission region, gamma-rays above 
10 ^^ as well as secondary electrons from vr^-decays effectively interact 
with the ambient photon and magnetic fields, and thus initiate IC 
and/or (depending on the strength of the B-field) synchrotron cascades. 
The cascade development stops when the typical energy of 7 rays is 
dragged to 1 TeV. Gamma-rays produced in this way can explain the 
observed TeV flux, if the acceleration power of 10^® eV protons is about 
IqS? ej-g/g_ xhe energy spectra of gamma-rays produced in this scenario 
are shown in Fig. 2. 

2.2.2. Proton-proton scenario 

Acceleration of protons to extremely high energies, E 10^® eV, is 
a key condition for effective photo-meson interactions. This model re¬ 
quires existence of strong magnetic field, B > 10^ G, in the compact 
region limited by a few gravitational radii. If the field close to the black 
hole is weaker, the interactions of protons with the ambient thermal gas 
become the main source of production of gamma-rays and electrons of 
“hadronic” origin. 

Protons can be accelerated to TeV energies also in the accretion 
disk, e.g. through strong shocks developed in the accretion flow. The 
efficiency of gamma-ray production in this case is determined by the 
ratio of accretion time R/v^ ~ 10 ^ — 10 *^ s (depending on the accretion 
regime) to the p-p cooling time, tpp ~ 1.5 x 10^(n/10® cm“^)“^ s. 
For any reasonable assumption concerning the density of the ambient 
thermal plasma and the accretion regime, the acceleration power of 
high energy protons should exceed Lp ~ 10 ^® erg/s in order to provide 
detectable fluxes of TeV gamma-rays (see Fig. 3). 

2.3. Curvature Radiation - Inverse Compton (CRIC) model 

The models of gamma-ray emission associated with accelerated protons 
provide rather modest efficiencies of conversion of the energy of acceler- 


AhaNer.tex; 2/02/2008; 16:58; p.4 



TeV radiation from GC 


5 



Figure 2. Broad-band Spectral Energy Distribution (SED) of gamma-rays (solid 
lines), neutrons (dash line) and neutrinos (dots) from Sgr A* due to interactions 
of ultra high energy protons with ambient photon and magnetic fields. Protons 
accelerated to energies eV in the regular magnetic field close to the gravitational 
radius R ^ Rg, propagate through the infrared emission region of size R = 107?g. 
The calculations correspond to two assumptions for the strength of the magnetic 
field in the region of the infrared emission: B — 0.1 G (a) and B = 10 G (b). The 
experimental points in X-ray and gamma-ray bands are from the compilation of 
Aharonian and Neronov (2005). 


ated protons to gamma-rays. The radiative energy loss rate of electrons 
is much higher, and therefore the models associated with accelerated 
electrons provide more economic ways of production of high energy 
gamma-rays. Obviously, these electrons should be accelerated to at 
least iiCmax ~ 10 TeV. This immediately constrains the strength of 
the chaotic component of the magnetic field; even under an extreme 
assumption that the acceleration proceeds at the maximum possible 
rate, {dE/dt)a_cc — ei?, one gets B < 10(iilmax/10TeV)“^ G. 

The requirement of particle acceleration at the maximum rate im¬ 
poses strong restrictions on the geometry of magnetic field and possible 
acceleration mechanisms. In this regard, acceleration in ordered electric 
and magnetic fields, e.g. by the rotation-induced electric, provides max¬ 
imum energy gain. Moreover, in the ordered field the energy dissipation 
of electrons is reduced to the curvature radiation loses; this allows 
acceleration of electrons up to £'e,max — 10^'^(i?/10G)^/^ eV. 

The curvature radiation peaks at Ccurv — 2 x 10®(iiie/10^^eV)^eV. 
The Compton scattering of same electrons leads to the second peak 
at much higher energies, ^ Eg ^ 10^^ eV (because the scattering 
proceeds in the Klein-Nishina limit). However, because of interactions 
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Figure 3. SED of the broad-band electromagnetic radiation initiated by p — p in¬ 
teractions in the accretion disk (from Aharonian and Neronov 2005). It is assumed 
that the accelerated protons with spectrum E~^ exp {—E/Ec) are injected into the 
thermal plasma of density 10® cm“®, and together with the accretion flow cross the 
region of the size R « 10i?g cm and fall under the black hole horizon after 10^ s. The 
following parameters have been assumed: (1) heavy solid curve: T = 2, high energy 
exponential cut-off at Ec = 100 TeV, total acceleration rate Lp = 5 x 10®® erg/s; 
(2) thin solid curve: F = 1, Ec = 1 TeV, Lp = 10"^® erg/s; (3) dashed curve: narrow 
(r = 0, Ec = 5 GeV) distribution of protons, Lp = 10"^° erg/s. For all three cases 
the magnetic field is assumed to be B = 10 G. 

with infrared photons, gamma-rays of energy exceeding 10 TeV can not 
freely escape the source. Synchrotron radiation and Compton scattering 
of the secondary (pair-produced) electrons lead to re-distribution of the 
initial gamma-ray spectrum. 

We call this scenario of production of Chrvature iJadiation and IC 
photons by electrons accelerated in regular magnetic/electric fields as 
CRIC model. Quantitative calculations of high energy radiation within 
framework of this model require a “self-consistent” approach which 
should take into account the radiation reaction force. An example of 
such self-consistent computation is shown in Fig. 4. 


3. Diffuse gamma radiation of runaway protons 

A significant fraction of protons accelerated near the black hole may 
escape the source and enter the surrounding dense gas environment. 
The interactions of these runaway protons with the interstellar medium 
lead to production of gamma-rays the luminosity of which could exceed 
the gamma-ray luminosity of the central engine. 
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Figure 4- Broad-band SED of radiation produced by electrons within the CRIC 
model (from Aharonian and Neronov 2005). Thin solid curve - gamma-ray pro¬ 
duction spectrum formed as superposition of the Curvature and inverse Compton 
emission components which accompany electron acceleration by the rotation-induced 
electric field within R = 2R^; heavy solid curve - the spectrum of gamma-rays mod¬ 
ified after the passage through the infrared source of size R = 10i?g. The strength 
of the regular magnetic field in the electron acceleration region is assumed B — 10 
G. The strength of the random magnetic field in the region of infrared emission is 
assumed B = 30 G. 


The flux of this radiation component depends not only on the 
density of the ambient gas and the injection rate of protons, but also 
on the speed of their propagation in the interstellar medium. The 
TeV radiation detected by HESS is localised within several angular 
minutes. This implies that, for the distance to the Galactic Center 
d ~ 8 kpc, the linear size of production region of gamma-rays can be 
as large as 10 pc. The density of the diffuse interstellar gas in this 
region is very high, n ~ 10^ cm“^. The lifetime of protons in this 
dense environment against nuclear p-p interactions is tpp ~ 5 x 10^ 
yr. Thus, in the case of absolute conhnement of protons, the gamma- 
ray luminosity after r\j years of continuous injection of protons 
achieves its maximum (“saturated”) level, = r^lTp, with an effi¬ 
ciency r/ ~ 1/3 (the fraction of energy released in neutral vr-mesons). 
However, the confinement time in the 10 pc region is rather limited; 
even in the Bohm diffusion regime the escape time of protons from 
this region is tesc ~ B?/vgC ~ 3 x 10®(E/100 TeV)“^(B/100 /rG) yr. 
Thus for any realistic diffusion coefficient D{E), TeV protons leave the 
region before they loose their energy in p-p interactions. This not only 
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reduces the gamma ray production efficiency but also, in the case of 
energy-dependent diffusion, modifies the energy distribution of protons 
established within the 10 pc region. 

The impact of the energy-dependent diffusion on the resulting 
radiation spectra is demonstrated in Fig.5. It is assumed that during 
10 ® years protons are injected (quasi)continuously into the interstellar 
medium of density n = 10® cm“®. The initial spectrum of protons is 
assumed in the form of power-law with an exponential cutoff, Q{E) = 
exp (—Fl/Flo)- The cutoff energy is fixed at Eq = 10^® eV, 
which is an obvious condition for effective production of gamma-rays 
to at least 10 TeV. The choice of the power-law index depends on 
the assumed diffusion coefficient, if one intends to explain the energy 
spectrum of gamma-rays detected by HESS, J{E) = (2.5 ± 0.21) x 
10“^^E“'" ph/cm^s TeV with T = 2.21 ±0.09 (Aharonian et al. 2004). 

The diffusion coefficient is assumed in the following form: D{E) = 
10^®(E/1 GeV)^K cm^/s. The values of k 1 and /3 0.5 correspond 

to the CR diffusion in the galactic disk. Of course, in the central region 
of the Galaxy one may expect significant deviation from the character 
of particle diffusion in “ordinary” parts of the galactic disk. In the case 
of effective confinement of protons, e.g. with n = 10“^ and (5 = 0.5, the 
escape time of multi-TeV protons is comparable with the characteristic 
time of p-p interactions. This prevents strong modification of the initial 
proton spectrum. Therefore, the injection spectrum of protons with 
power-law index a = 2.2 fits quite well the observed TeV spectrum 
(curve 2 of Fig.5). Because of the effective confinement, the required 
injection rate of protons is rather modest, lEp = 7 x 10®® erg/s. 

In the case of faster diffusion, the spectral index of injection is 
determined as a ~ T — /3. For example, for n = 0.15 and (5 = 0.3 
(a diffusion regime which corresponds to the Kolmogorov type turbu¬ 
lence) , the 7 -ray observations are well explained assuming the following 
parameters a = 1.9 and IFp = 7.5 x 10®'^ erg/s (curve 1 in Fig.5). 

Finally, one should note that the injection spectrum remains un¬ 
changed with very fast escape. For example, for n = 10“^ and (5=1 
(an energy-dependence close to the one expected in the Bohm diffusion, 
but with much larger diffusion coefficient), the particles propagate in 
the diffusion regime until energies of several TeV. At higher energies 
they escape the source almost rectilinearly on timescales R/c ~ 30 
years. This implies that while at low energies the protons inside the 10 
pc region suffer significant spectral deformation, a' ^ a ± /3, at very 
high energies the shape of the initial spectrum is essentially recovered. 
Such an interesting modulation of the proton spectrum is reflected in 
the resulting gamma-ray spectrum (curve 3 in Fig.5). In this case the 
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Figure 5. Energy spectra of diffuse gamma rays expected from the central 10 pc 
region of the Galactic Center within different assumptions concerning the diffusion 
coefficient and the proton injection spectrum (see the text). 


requirement to the injection power of protons is higher than in the 
previous cases, Wp = 10 ^® erg/s. 

The TT^-decay radiation of protons is always accompanied by syn¬ 
chrotron radiation of secondary electrons - the products of charged 
vr-mesons. In the case of extension of the proton spectrum to 10^® eV, 
and for the magnetic field exceeding 100 /tG, the spectral energy distri¬ 
bution (SED) of synchrotron radiation of secondary electrons peaks in 
the X-ray domain. This seems an attractive mechanism for explanation 
of the diffuse X-ray emission of the Galactic Genter, given the seri¬ 
ous problems of interpretation of this radiation within the “standard” 
(thermal and nonthermal) models (Muno et ah, 2004). However, for a 
relatively flat SED of TeV gamma-rays (like the one observed by HESS) 
the energy flux of X-rays is always less, by a factor of 3 to 10, than the 
energy flux of 7 -rays. Therefore for the spectra shown in Fig. 5 only 
a small fraction (10 per cent or so) of the observed X-ray flux can be 
contributed by secondary electrons. 

Formally, the flux of the secondary-electron-synchrotron compo¬ 
nent can be significantly increased, without getting in conflict with 
the TeV data, if one assumes very hard spectrum of protons, e.g. with 
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Figure 6. The broad-band spectrum of radiation initiated by interactions of ul¬ 
tra-high energy protons in the central 10 pc region. The injection spectrum is 
assumed in the form of power-law with a = 1.5 and low and high energy cutoffs at 
= 2 X lO^'* eV and E\ = 1.5 x 10^® eV. The proton injection rate is assumed 
b = 1.5 X 10®® erg/s to match the reported X-ray (Chandra - Muno et al. 2004, 
INTEGRAL - Belanger et al. 2004) and TeV fluxes (Aharonian et ah, 2004). The 
diffusion coefficient is characterised by /3 = 0.3 and k = 10“®. The number density 
of the gas and magnetic field in the central 10 pc region are assumed n = 10® cm“® 
and B = 1 mG, respectively. The effect of spectral modulation due interactions with 
the 2.7 K CMBR is not shown. Note that for the distance to the GG the absorption 
becomes noticeable (~ ~ 1/3) only at energies around 10^® eV. 


power-law spectral index a < 1.5. Such a spectrum cannot explain the 
TeV data. However, assuming that the spectrum of protons extends to 
Eo r\j 10^® eV, a new TeV gamma-component, due to the synchrotron 
radiation of the secondary ultrahigh energy electrons, could dominate 
over the vr’^-decay component. The possibility of producing high fluxes 
of X-rays and gamma rays through the synchrotron radiation of sec¬ 
ondary electrons is demonstrated in Fig. 6. Note that although the 
parent protons interact with the ambient gas throughout the central 
10 pc region of the Galactic Center (thus the radiation can be dubbed 
as diffuse), the observer would detect a point-source like TeV signal 
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centered on Sgr A*. The reason is that the ultra high energy protons 
injected into the surrounding medium propagate radially (like pho¬ 
tons) without significant deflection in the interstellar magnetic field. 
Since the characteristic energy of parent protons, which are respon¬ 
sible for the secondary synchrotron X-rays, is smaller by a factor of 
(TeV/keV)^/^ ~ 10^ — 10^, the angular size of X-rays should be larger. 
This scenario predicts a tendency of decrease of the angular size of the 
X-ray signal with increase of the photon energy. This hypothesis can 
be inspected also by detection of direct gamma-rays from photome¬ 
son interactions up to ~ 10^^ eV. Above several TeV this radiation 
component dominates over the secondary synchrotron radiation. 


4. Summary 

The origin of TeV radiation detected from the direction of the Galactic 
Center is not yet established. The 3-arcmin upper limit on the source 
size, assuming Gaussian distribution of the source brightness (Aharo- 
nian et al. 2004), implies that several objects, in particular the Dark 
Matter Halo, the young SNR Sgr A East, the central compact source 
Sgr A*, as well as the entire central diffuse region filled by dense molecu¬ 
lar clouds and cosmic rays, are likely candidates for TeV emission within 
the central 10 pc region. In fact, each of these sources may contribute 
significantly into the observed gamma-ray flux. In this paper we present 
the results of our study of several gamma-ray production scenarios with 
both direct and indirect links to the massive black hole in the Galactic 
Center, Sgr A*. 

The results shown in Figs. 2-4 demonstrate that at least three 
gamma-radiation scenarios, which take place in vicinity of the massive 
black hole, can explain the TeV observations without apparent conflicts 
with observations of Sgr A* at lower frequencies. The hadronic mod¬ 
els based on photo-meson or pp interactions cannot provide efficiency 
higher than 0.1 per cent, therefore they require an acceleration rate 
of protons of about 10^® erg/s or larger. Although this is larger than 
the total electromagnetic luminosity of Sgr A*, it is still acceptable for 
a black hole of mass 3 x 10 ®Mq. On the other hand, the electronic 
model CRIC, based on the curvature and synchrotron radiation chan¬ 
nels, allows an economic way of conversion of energy of electrons to 
gamma-rays. 

The pp hadronic model predicts X-ray fluxes comparable with X- 
ray observations in high state, but an order of magnitude higher that 
in the quiescent state. The photomeson and CRIC models predict an 
order of magnitude lower X-ray fluxes. On the other hand, the photome- 
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son model predicts detectable fluxes of ultrahigh energy gamma-rays 
and neutrons. The CRIC model predicts high gamma-ray fluxes at 
MeV/GeV energies which should be easily detected by GLAST. 

The gamma-radiation in all three models is generally expected to 
be variable on timescales as short as 1 hour. Therefore the detection of 
a variable component of radiation on such short timescales would be a 
strong argument in favour of gamma-ray production near the massive 
black hole. On the other, although the lack of variations of the TeV flux 
could be naturally interpreted as gamma-ray production in extended 
regions, it cannot be used as a decisive argument against the black-hole 
origin of TeV emission. 

In addition to the gamma-rays emitted in compact regions in 
the vicinity of the massive black hole, one should expect also a dif¬ 
fuse (extended) component of radiation associated with interactions 
of the runaway protons with the surrounding dense interstellar gas. 
The relative contribution of this component to the total TeV flux, as 
well as the spectral and angular characteristics of vr^-decay gamma- 
rays significantly depend on the rate of injection of protons by Sgr 
A* into the interstellar medium, as well as on the regime of (energy- 
dependent) diffusion of protons in the dense central region. For certain 
combinations of principal model parameters the diffuse gamma-rays 
from pp interactions can satisfactorily explain both the absolute flux 
and energy spectrum of TeV radiation reported by HESS. Although this 
component can be extended well beyond the Galactic Center region, 
because of enhanced gas density in the central 10 pc region one should 
expect a bright gamma-ray core within several arcminutes around Sgr 
A*. In the case of extension of the spectrum of protons to E ~ 10^® eV, 
one may expect another diffuse component related to the synchrotron 
radiation of ultrahigh energy secondary electrons (from yr^-decays). 
Since the gyroradius of parent protons is comparable or exceeds 10 pc, 
the observer would detect this TeV component as a point-like source. 
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